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Introduction 
Magneto-fluid  dynamics  and  plasma  physics  are  a 
continuing  source  of  nonstandard  and  unusual  mathematical 
structures  (e.g.  see  Ref s. [1] , [2] ,  [5] ) .   Once  such 
question  arises  in  a  hyperbolic  system  in  which  the  time- 
like direction  is  periodic  or  re-entrant.   The  problem 
which  we  consider  here  is  that  of  a  composite  system, 
partly  elliptic  and  partly  degenerate  hyperbolic  (i.e. 

with  real  characteristics  that  are  lines  instead  of  cones),*  in 
a  toroidal  domain „   It  has  been  found  in  this  case  that  it 

is  impossible  to  formulate  a  general  boundary  value  problem 

with  smooth  solutions  [4].   There  is  an  open  question  of 

the  possibility  of  a  formulation  in  terms  of  appropriate 

weak  solutions;  also  an  Important  practical  question  is 

the  Interpretation  of  the  large  literature  of  formal 

calculations  which  give  explicit  answers  where  none  exist  I 

The  Differential  Equations 
The  problem  is  the  static  balance  of  a  perfectly 
conducting  plasma  in  equilibrium  with  a  magnetic  field. 
In  a  macroscopic,  scalar  pressure,  MHD  version  we  have 

Vp  =  J  X  B 

(1) 

J  =  curl  B,  div  B  =  0 


The  characteristics  of  this  system  are  one  elliptic 
cone  (as  in  the  special  case,  curl  B  =  0)  plus  the 
magnetic  lines  counted  twice  as  degenerate  real 
characteristics  (cof.  streamlines  in  fluid  flow). 

It  is  frequently  convenient  to  parameterize 
div  B  =  0  locally  in  terms  of  two  stream  functions, 

B  =  Va  X  vp  (2) 

One  of  the  two  sets  of  real  characteristics  can  some- 
times be  eliminated  by  considering  p  =  p{ot,  P)  to  be  a 
known  function;  i.e.  by  assigning  a  value  of  p  to  each 
magnetic  line--but  this  may  be  troublesome  in  the  large. 

In  an  anisotropic,  guiding  center  version,  instead 
of  p  =  p(ct,  P)  we  introduce  p  =  p(a,  P,B)  where  the  absolute 
value  of  the  field  strength,  B,  is  taken  as  the  third 
coordinate  instead  of  arc  length  along  a  magnetic  line. 
We  obtain  the  system  [3] 


Is  '-  +  3^P  'f'  =  K  X 


B 


K  -  curl  (6B),  div  B  -  0  (3) 

^  -  ^        B  ^B 
If  p(a, p^B)  is  considered  to  be  given,  the  family  of 
magnetic  field  lines  remains  as  one  real  characteristic, 
and  there  is  an  additional  nondegenerate  cone  which  is 
elliptic  provided  that  6  >   0   and  S((5B)/BB  >  0, 

bTb   -^  (2|) 

^  ,1   hp^         , 
^^B  m'         ^ 


If  these  inequalities  (which  refer  to  the  'firehose'  and 
'mirror'  instabilities  respectively)  are  violated, (3)  does  not 
give  a  well-p'jsed  boundary  value  problem  for  B  [5]. 

In  terms  of  conventional  physical  variables,  the 
quantity  we  call  p  is  p^^,  the  parallel  component  of  the 
anisotropic  pressure  tensor;  p  =  -B  S(p/B)/c)B  is  the 
perpendicular  component,  and  6-1  =  (p  -p^  )/B  is  a  measure 
of  the  anlsotropy.   The  physical  condition  that  p  be 

positive,  S(p/B)/^B  <  0,  implies  the  firehose  condition 

2 
when  p  <  B  and  is  implied  by  it  otherwise. 

With  axial  symmetry,  writing 

B  -  |x  V-^  +  uBq  (5) 

we  find  that  both  real  characteristics  integrate  out 
explicitly,  one  through  p  =  p(^, B)  and  the  other  through 

I  er^B^  -  g(^)  (6) 

to  give  a  nominally  elliptic  equation 

^  +  r  J-i^  1^-)  +  ^(1^  4^-  +  ^  ^) 
:,  2  ^  ^  ^^r  B^r^  dMz  ^z  ^  Br' 
oz 


-^^g'(^)^-#if  -  ° 


(7) 


[the  given  function,  6    =  6(^,B),  where  B  =  \^f\^/r     -h   2g/6r  , 
determines  d  implicitly  in  terms  of  r,^,  and  |V^|]. 

For  a  more  detailed  discussion  of  the  guiding-center 
equations,  see  Refs.[33  ai^d  [5]. 


Boundary  Value  Problems 

A  plausible  theory  can  be  given,  based  In  part  on 
proved  existence  theorems,  for  a  tubular  domain  with  end 
surfaces  on  which  B^  >  0  and  B^  <  0  respectively  and  a 
lateral  surface  where  B  =  0.   The  domain  is  assumed  to 
be  simply  covered  by  the  direction  field  B.   Corresponding 
to  the  elliptic  characteristic  cone  we  specify  the  normal 
component,  B  ,  subject  to  pB  dS  =  0,  on  the  entire  boundary. 
Corresponding  to  the  real  characteristics,  we  specify 
two  pieces  of  initial  data  at  one  end  of  the  tube  (or  one 
at  each  end);  e.g.  we  may  specify  p  and  J^.   In  addition 
to  qualitative  arguments  such  as  counting  characteristics, 
we  confirm  a  consistent  variational  description  with  im- 
posed or  natural  boundary  data,  we  verify  actual  existence 
in  the  case  of  axial  symmetry,  and  we  are  able  to  construct 
plausible  iteration  procedures  which  rationally  account 
for  all  boundary  data  and  which  can  be  shown  to  converge 
with  suitable  restrictions  in  axial  symmetry  (and  should 
probably  converge  more  generally);  see  Ref.[6]<. 

Similar  conjectures  have  been  made  for  a  toroidal 

domain  [6]  (and  without  even  recognizing  them  as  conjectures 

[7]),assuming  that  the  domain  is  covered  by  a  family  of 

nested  toroidal  pressure  surfaces.   In  the  absence  of  an 

'initial'  surface  which  intersects  each  magnetic  line, 

two  pieces  of  information,  e.g.  p(V')  and  the  rotation 

number  cL{f),    are  assigned  to  each  flux  surface.   The 

plausibility  arguments  such  as  characteristic  counting, 

variational  analysis,  and  actual  proved  examples  in  axial 

syinmeti'y  are  equally  valid  for  the  torus.   In  addition, 
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the  fact  that  the  tangent  field  B  satisfies  an  elliptic 
equation  on  each  pressure  surface  refines  the  counting 
by  demanding  that  two  periods  be  assigned  to  each  pressure 
surface o 

A   small  shadoiv  was  cast  even  at  the  beginning,  by 
the  fact  that  iterations  that  were  plausible  in  the  open 
tube  could  not  be  made  to  work  for  the  torus.   After  a 
series  of  disturbing  irregularities  in  expansions  (e.g. 
see  [8]),  it  was  finally  recognized  that  axial  (or  other) 
symmetry  is  not  just  a  convenience  for  calculation,  but 
it  is  necessary  for  existence  [4]. 

The  situation  is  qualitatively  similar  to  the  question 

of  existence  of  flux  surfaces  for  harmonic  vector  fields, 

curl  B  =  0,  in  general  toroidal  domains,  but  it  is  not 

the  same.   The  possibility  of  taking  J  ^  0  allows  more 

flexibility.   For  example,  given  an  arbitrary  set  of 

nested  flux  surfaces,  the  necessary  condition,  Vp.B  x  curl  B  =  0 

div  B  =  0,  can  be  satisfied  together  with  arbitrarily  given 

p(^)  and  a(^).   The  problem  of  finding  equilibria  is,  also 

connected  with  (but  again  is  not  a  simple  consequence  of) 

the  pathology  of  particle  orbits  in  a  toroidal  domain. 
It  is  most  directly  related  to  the  very  singular  guiding 

center  limit  which  permits  infinite  currents  to  develop  [h]. 
Two  problems  immediately  present  themselves.   The  first 

is  how  to  alter  the  mathematical  formulation  in  order  to 

restore  existence.   The  second  is  how  to  interpret  the  large 

literature  of  formal  calculations  in  specific  geometries 

which  give  explicit  results-"but  even  in  cases  where  no 

correct  answer  exists. 


With  regard  to  restoring  existence,  the  microscopic 
model  which  leads  to  the  anisotropic  system  (5)  Is  help- 
ful In  associating  mathematical  nonexistence  with  actual 
physical  pathology  [4]  (In  some  cases,  this  has  been  ob- 
served experimentally).   In  other  words,  there  are  cases 
In  which  we  cannot  expect  help  from  any  Improved  physical 
model;  nonexistence  of  a  stationary  solution  Is  the  correct 
physical  predlctlono   Qualitatively,  we  are  able  to  draw 
the  conclusion  that  In  the  limit  of  small  Larmor  radius,  a 
smooth  pressure  profile  must  be  associated  with  pathological 
current  (J  Is  nonlntegrable ).  On  the  other  hand,  assuming  that 
J  remains  bounded  Implies  that  Vp  is  pathologlcal-~at 

best  measur.s'able,  but  never  continuous.   In  other  words, 
p  and  B  cannot  both  be  smooth.   We  are  presented  with  an 
unusual  situation  In  partial  differential  equations  where 
weak  solutions  are  not  adopted  as  a  mathematical  convenience  but 
at  the  Insistence  of  the  equations. 

It  Is  Interesting  to  note  that  the  theory  of  Arnold  [9] 
In  this  problem,  although  technically  correct.  Is  vacuous 
because  It  assumes  analytlclty  of  p  and  of  B. 

The  question  of  the  existence  of  weak  solutions  to 
the  systems  (1)  or  (3)  Is  open  (although  some  preliminary 
analysis  has  been  done  In  weak  perturbations  of  a  given 
strong  solution  P-O]).   To  find  a  weak  solution,  it  is 
first  necessary  to  specify  data  such  as  p(^)  in  situations 
where  the  variable  f   is  not  well-defined.   We  may  consider, 
in  accordance  with  the  theories  of  Arnold  [u]  and  Moser  [1^, 
that  there  is  a  family  of  surfaces,  tp  =   const.,  of  positive 
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measure,  together  with  gaps  (ergodlc,  islands,  etc. )  surrounding 
regions  of  rational  rotation  number,  a.   We  can  require 

that  p  be  constant  in  each  gap.   It  Is  also  plausible 

that  a  is  constant  throughout  each  gap  (it  is,  of  course, 

constant  in  any  Island).  A   smooth  function  p(a)  might 
therefore  be  specified  as  a  constraint  which  could  be 
well-defined  even  for  weak  solutions.   Such  an  analysis 
has  not  yet  been  done. 

The  second  problem  mentioned  above  is  the  interpreta- 
tion of  formal  calculations,  usually  expansions,  which 
present  ostensibly  valid  equilibria,  complete  with  flux 
surfaces.   In  some  cases  the  pathology  is  hidden  in  a 
higher  order  of  expansion  (e.g.  if  a  gap  has  thickness  e  ). 
In  other  cases,  the  pathology  appears  to  be  transcendental 
and  is  not  seen  to  any  order. 

The  salvation  of  such  calculations  rests  on  the  fact 
that  the  time-dependent  MHD  initial  value  problem  is  well- 
posedo   Pathology  enters  only  when  we  ask  the  improper 
question  of  the  existence  of  a  static  (or  even  stationary) 
solution. 

The  existence  of  static  solutions  is  a  very  great 
convenience  in  classical  analysis.   For  example,  almost 
all  stability  analysis  introduces  a  small  perturbation 
from  a  postulated  equilibrium.   We  can  justify  such  MHD 
stability  calculations  as  follows. 

The  equations  of  motion  are  hyperbolic  with  character- 
istics determined  by  the  Instantaneous  state,  viz.  p(x),  u(x), 
and  B(x).   Suppose  we  are  given  truncated  power  series. 


p(^)  =  E  e^p   ,  B^"^  -  2  e^B   ,  such  that 
o    ""^  o    ■"" 

Vp(")  ^  J^^^  X  B^^^  +  e^^^R^+i  (8) 

and  R   (x, e)  Is  a  knovm  smooth  function  of  x.   The  problem 
is  to  confirm  whether  the  solution  of  the  hyperbolic  system 
with  initial  values  u  =0,  p  =  p   ,  B  =  B^  Remains  near 
the  initial  state  for  an  appreciable  time. 

As  a  simplification,  consider  the  symmetric  hyperbolic 
system  obtained  by  a  conventional  linearization  of  the 
equation  of  motion  about  the  state  u  =  0,  p=p^   ,B  =  B^  '. 
This  system  is  inhomogeneous  through  R  -,,  has  zero  initial 
values,  and  has  fixed  characteristics  determined  by  p 
and  B^^  .   The  solutions  will  remain  bounded  (in  a  suitable 
norm)  if  the  system  is  stable  in  a  conventional  sense.   The 
crucial  point  is  that  the  theory  of  the  hyperbolic  system 
is  standard  when  linearized  about  p^^  and  B^^^  whether 
or  not  this  state  corresponds  to  a  true  static  equilibrium. 

At  present,  the  information  available  is  more 
limited.   The  stability  analysis  is  itself  carried  out 
as  an  expansion  in  e.   Since  in  the  m.ost  common  case 
there  is  a  highly  degenerate  eigenvalue,  ij  =  0,  which  is 
spread  by  the  parameter  e,  even  the  linear  program. 
mentioned  above  has  not  been  carried  out„   But  it  is  at 
least  possible  to  point  to  a  conventional  mathematical 
framework  within  which  we  could  determine  the  significance 
of  formal  expansions  of  nonexistent  equilibria. 

8 


There  Is  one  sense  in  which  the  formal  expansions 
do  reflect  a  shadow  of  the  nonexistence  of  a  legitimate 
solution.   In  common  with  many  asymptotic  expansions,  the 
determination  of  a  given  order  of  expansion  is  not  com- 
plete until  partial  information  is  extracted  from  a 

higher  order.    In  the  present  case,  we  find  that 

n4-2 
field  perturbations  of  order  e    and  sometimes  even 

n+4 
e    react  back  on  the  equilibrium  configuration  to 

order  e  ,   The  (nonexistent)  calculated  configuration 

is  therefore  very  sensitive  to  perturbations  [IJ)]. 


The  Free  Boundary 
There  is  one  special  case  in  which  the  nature  of  the 
weak  equilibrium  can  be  ascertained.   The  MHD  free  boundary 
is   a  flux  surface  separating  a  harmonic  field,  curl  B  =0, 
div  B  =  0,  from  a  field-free  region  at  constant  pressure  p,  and 
is  mathematically  identical  to  the  classical,  fluid  free 

boundary  in  virtue  of  the  boundary  condition,  p  =  (l/2)B  or 

2 

B  =  const. 

There  is  a  large  literature  concerning  two-dimensional 
and  axially  symmetric  free  boundaries.   One  attack  which 
we  can  hope  to  extend  to  three  dimensions  is  the  Inverse 
method.   Given  an  analytic  curve  on  which  f  -  0   and 
hf/^n   =  1,  we  employ  the  Cauchy-Kovalevsky  theorem  to 
solve  an  initial  value  problem  for  A-^  =  0  in  the  neigh- 
borhood of  the  curve. 

In  three  dimensions,  there  is  no  equation  equivalent 
to  A^  =  0.   To  use  the  potential,  B  =  V0,  A0=  0,  we  must 
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first  find,  on  a  given  analytic  surface,  a  scalar  0 
with  constant  surface  gradient,  |V0|  =1.   The  Cauchy 
Initial  values  will  then  be  0  and  B^/Sn  =  0. 

This  Is  an  essentially  hyperbolic  problem  on  the 
initial  surface.   The  Hamllton-Jacobl  wave  fronts, 
0  =  const.,  are  orthogonal  to  geodesic  rays  (magnetic 
field  lines  on  the  surface).   In  other  words,  if  on  a 
given  analytic  toroidal  surface,  we  are  given  an  analytic 
single-valued  field  of  geodesies,  we  can  then  construct 
the  orthogonal  trajectories  and  assign  a  multivalued 
function  0  with  single-valued  V0  and  |V0|  =1.   Existence 
of  a  three-dimensional  field  B  =  V0  in  the  neighborhood 
of  the  surface  then  follows. 

The  question  reduces  to  that  of  finding  a  geodesic 
field--thls,  of  course,  is  not  alvmys  possible.   If  we 
treat  the  problem  by  perturbation  about  an  elementary 
surface  (e.g.  a  cylinder  with  identified  ends  or  an 
axlally  symmetric  surface)  on  which  we  know  of  the  exist- 
ence of  a  geodesic  field  for  each  value  of  the  rotation 
number  o-,    the  problem  is  reduced  to  an  application  of 
the  twist   mapping  theorem  of  Moser  and  Arnold.   In 
other  words,  there  will  be  a  geodesic  field  on  any  sur- 
face which  is  close  to  an  'elementary'  surface  for  a 
set  of  a  of  positive  measure  which  avoids  the  neighbor- 
hood of  rational  a. 

Prom  this  result,  wa  can  conjecture  the  structure 
of  the  solution  of  the  direct  problem.   In  two  dimensions, 
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for  a  given  fixed  outer  boundary,  there  exists  a  one 
parameter  family  of  Interior  free  boundaries,  e.g.  para- 
metrized by  area.   In  three  dimensions,  let  us  choose  a 
as  the  parameter  Instead  of  area  (or  volume).   We  can 
expect  to  find  a  family  of  free  boundaries  corresponding 
to  a  set  of  positive  measure  of  a,  but  not  a  complete 
family,  for  each  value  of  the  area,  as  In  two  dimensions, 
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Magneto-fluid  dynamics  and  plasma  physics  are  a  continuing  source  of 
nonstandard  and  unusual  mathematical  structures.   Once  such  question 
arises  in  a  hyperbolic  system  in  which  the  time-like  direction  is 
periodic  or  re-entrant.   The  problem  which  we  consider  here  is  that 
of  a  composit  system,  partly  elliptic  and  partly  degenerate  hyper- 
bolic (i.e.  with  real  characteristics  that  are  lines  instead  of 
cones),  in  a  toroidal  domain.   It  has  been  found  in  this  case  that  it 
is  impossible  to  formulate  a  general  boundary  value  problem  with 
smooth  solutions.   There  is  an  open  question  of  the  possibility  of 
a  formulation  in  terms  of  appropriate  weak  solutions;  also  an 
important  practical  question  is  the  interpretation  of  the  large  litera- 
ture of  formal  calculations  which  give  explicit  answers  where  none 
exist  I 
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